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Abstract — We present two interference alignment tecli- 
niques sucli tliat an opportunistic point-to-point multiple 
input multiple output (MIMO) link can reuse, without 
generating any additional interference, the same frequency 
band of a similar pre-existing primary link. In this 
scenario, we exploit the fact that under power constraints, 
although each radio maximizes independently its rate by 
water-filling on their channel transfer matrix singular 
values, frequently, not all of them are used. Therefore, 
by aligning the interference of the opportunistic radio it 
is possible to transmit at a significant rate while insuring 
zero-interference on the pre-existing link. We propose a 
linear pre-coder for a perfect interference alignment and 
a power allocation scheme which maximizes the individual 
data rate of the secondary link. Our numerical results show 
that significant data rates are achieved even for a reduced 
number of antennas. 

I. Introduction 

We consider the case of radio devices attempting to 
opportunistically exploit the same frequency band being 
utilized by licensed networks under the constraint that 
no additional interference must be generated. This can be 
implemented for example by assuming that opportunistic 
users are cognitive radios |T1|, f2\. Typically, cognitive 
radios temporally exploit the unused frequency bands, 
named white-spaces, to transmit their data. This clearly 
improves the spectral efficiency since more users are 
allowed to co-exist in the same bandwidth. However, a 
higher spectrum efficiency could be attained by simul- 
taneously allowing the opportunistic radios to transmit 
with the licensees if no harmful interference is gener- 
ated. In this case, interference alignment (lA) has been 
identified as a powerful tool to achieve such a goal. 
This technique, was initially introduced in |f3l, HI. lA 
allows a given transmitter to partially or completely 
"align" its interference with unused dimensions of the 
primary terminals [5|. The concept of dimension could 



Figure 1. Two-user MIMO interference channel. 



be associated with a specific spatial direction, frequency 
carrier or time slot H, ||6l. An extensive study has 
been conducted H, Q, Q to estimate the number of 
interference-free dimensions a given radio might find to 
transmit when several radio systems co-exist. In [4J and 
[6] several lA schemes to exploit such dimensions are 
proposed. Similarly, in fSl, a linear pre-coder based on 
Vandermonde matrices allows an orthogonal frequency 
division multiplexing (OFDM) radio to co-exist with 
similar pre-existing terminals without generating any 
additional interference. The idea is to exploit the re- 
dundancy of the OFDM cyclic prefix and frequency 
selectivity of the channel. 

In this study, we propose a novel interference align- 
ment technique for the secondary users exploiting the 
fact that under a power-limitation, a primary user which 
maximizes its own rate by water-filling on its MIMO 
channel singular values, might leave some of them 
unused, i.e. no transmission takes place along the corre- 
sponding spatial directions. These unused directions may 
be opportunistically utilized by a secondary transmitter. 



since its signal would not interfere with the signal sent 
by the primary transmitter. We present a linear pre-coder 
which perfectly aligns the interference generated by the 
secondary transmitter with such unused spatial direc- 
tions. Similarly, we present a power allocation scheme 
based on the water-filling idea which maximizes the 
individual data rate of the opportunistic radio. Simulation 
results show that a significant data rate can be achieved 
by the secondary link following our approach. 

II. System Model 

Notation: In the following, matrices and vectors are 
denoted by boldface upper case symbols and boldface 
lower case symbols, respectively. The i^^ entry of the 
vector X is denoted x{i). The entry corresponding to the 
i^^ row and j*^ column of the matrix X is denoted by 
X{i,j). The A^-dimension identity and null matrix are 
represented hy and Oat, respectively. The Hermitian 
transpose is denoted (•)^, and the expected value is 
represented by the operator E [.]. 

We consider two point-to-point unidirectional links 
simultaneously operating in the same frequency band and 
producing mutual interference as shown in Fig. [T] 
Both transmitters are equipped with Nt antennas while 
both receivers use A^^^. antennas. The first transmitter- 
receiver pair, i.e. Txi and Rxi, is a primary link licensed 
to exclusively exploit a given frequency band. The pair, 
Tx2 and Rx2 is an opportunistic link exploiting the same 
frequency band subject to the constraint that no addi- 
tional interference must be generated over the primary 
system. Note that no cooperation between terminals 
is allowed, i.e. transmitters do not share or exchange 
any signal before transmitting. Therefore, the multiple 
access interference (MAI) is considered as additive white 
Gaussian noise (AWGN). 

The channel transfer matrix from transmitter j G 
{1, 2} to receiver i G {1, 2} is denoted Hij G C^-^^', 
where the entries of Hij are independent and identically 
distributed (i.i.d) complex Gaussian circularly symmetric 
random variables. The channel matrices are supposed 
to be fixed for the whole transmission duration. This 
correspond to assuming (static) Gaussian links. But 
our analysis readily extends to the case of slow-fading 
channels by assuming the channels to be constant over 
each data block. Regarding the channel state informa- 
tion (CSI) conditions, we assume the primary terminals 
(transmitter and receiver) to only have perfect knowledge 
of the matrix Hi i. On the other hand, the secondary 
terminals have perfect knowledge of all the channel 
transfer matrices Hij, for every i and j G {1,2}. 
Although unrealistic, this condition provides us with an 
upper bound on the achievable rate of the secondary user. 
It can however be met in practice in the TDD (Time 
Division Duplex) mode if the secondary user exploits 



opportunistically the training sequences and signaling 
communication between the primary devices. 

Following a matrix notation, the primary and sec- 
ondary received signals can be written as 

f yi\_f Hi^2 \ ( Visi 

[y^J I H2,i H2,2 J I V2S2 / V ^2 r 

(1) 

where the vectors Si G C^' ^ ^ and rii G C^'' ^ ^ represent 
the transmitted symbols and an AWGN process with zero 
mean and covariance matrix u^/at^ for the i*^ link. For 
all z G {1,2} the matrices Vi G C^'^^' represent 
the linear pre-coders used for interference alignment. 
Furthermore, at each receiver, the input signal is linearly 
processed with the matrix Fi G C^''^^''. The signal 
at the output of the linear filter i is rj = Fiy^. Both 
matrices Vi and Fi are described later on. The power 
allocation matrices are defined as the input covariance 
matrices P, = E [s^sf ] G for the i*^ 

transmitter. The power constraints are 

V* G {1, 2} , Trace (F^P^Ff ) ^ Pi,max, (2) 

where Pi,max is the maximum transmit power level for 
the i^'^ transmitter. Without loss of generality, we assume 
identical maximum transmit powers for all the terminals 

i.S. Vl G {1, 2}, Pi^max — Pmax- 

III. Interference Alignment Strategy 

In this section we focus on the study of the pre-coding 
Vi and post-processing Fi matrices. Suppose that the 
primary terminals completely ignore the presence of the 
opportunistic transmitter. Hence, in order to maximize its 
own data rate, the primary transmitter follows a water- 
filling power allocation as in the single-user case lilOJ . 

A. Primary link design 

Under the assumption that the channel matrix -ff 1,1 
is known at the receiver and transmitter, the primary 
terminal chooses its pre-coding Vi and post-processing 
Fi matrices in such a way that their channel transfer 
matrix is diagonalized, i.e. Vi and Fi = Ui satisfy 
the singular value decomposition Hi 1 = UiAiV^ , 
where Ui G C^-^^- and Vi G C^'^^* are unitary 
matrices and Ai G (M+)^''^^* is a diagonal matrix 
which contains m.m{Nr,Nt) non-zero singular values, 
Ai, . . . , Ajnin(Af,,Aft)- Thus, the received signal after linear 
processing, ri, can be written as 

ri = Fiyi = Aisi + n[, (3) 

where n[ = U^ni is an AWGN process with zero 
mean and covariance matrix (j'^In^ - Then, the achievable 
rate of the primary user is maximized by the power 



allocation matrix Pi which is a solution to the following 
optimization problem 

maximize log2 |/Ar, + ^-Ffi^iFiPiVf i3"f^i| 
subject to Trace (Pi) ^ Pmax- 

The solution to (|4]| is the classical water-filling algorithm 
lITOl . Following this approach, the optimal power alloca- 
tion matrix is a diagonal matrix with entries 



yn£{l,...,Nt}, Pi(n,n) 



(5) 



with [p]~^ = max {0,p). The constant /3 is a Lagrangian 
multiplier that is determined to satisfy 



N 

i=i 



Pn 



B. Secondary link design 

Depending on the channel singular values 
Ai, . . . , Aniin(7v^ TVt)' t^c powcr allocation matrix 
Pi might contain zeros in its main diagonal. A zero 
power allocation for a given singular value means that 
no transmission takes place along the corresponding 
spatial direction. This means that the secondary terminal 
can align its transmitted signal with the unused singular 
modes such that it does not interfere with the signal 
transmitted by the primary user. If one converts the 
spatial problem into the frequency one, the result is 
similar to the cognitive scenario where the secondary 
would opportunistically use the unexploited frequency 
modes. The main difference here lies in the fact that in 
the spatial domain, there is no universal precoder which 
diagonalizes the basis of all the devices whereas this is 
the case in the frequency domain with the use of the 
FFT 

As a consequence, in the spatial domain, the cor- 
responding orthogonality condition (such that the sec- 
ondary user generates no interference on the primary 
link) is given by 



U?H,oV 



aPi, 



(6) 



where the matrix Pi is a diagonal matrix with entries 

,2 



VnG {!,..., iVt}, Pi(n,n) 



(7) 



such that the condition Pi Pi = Oat^ always holds. It 
can be easily verified since both matrices are diagonal. 
Additionally, the constant a is chosen to satisfy the 
power constraints Q with i = 2. 

Assuming that perfect estimates of Hi i and H12 
are available at the secondary transmitter, the secondary 
precoder (when the inverse of Hi^2 exists) is given by: 



(8) 



For the case where > Nt, i.e. the receiver has 
more antennas than the transmitter, it is still possible to 
obtain the pre-coding matrix by using the Moore-Penrose 
pseudo-inverse of Hi ^2, 



a 



[Hi2Hl^2) H12U1P1 



(9) 



Once the pre-decoder V2 has been adapted to satisfy 
dUl or (|9l) at the secondary transmitter, no additional 
interference impairs the primary user. However, the 
secondary receiver still undergoes the interference from 
the primary transmitter. Typically, this effect is a colored 
noise with covariance 

Q e £N,xN, due to the channel 
H2,i and the pre-coder Vi. Here, 



Q 



H2,iViPiVi H2^i 



N,.- 



(10) 



Hence, the received signal y2 can be whitened by using 
the matrix F2 = Q^^ , to obtain r2 = P2Z/2' such that 



r2 = Q ^2H2,2V2S2 + ^2, 



(11) 



where n'g = Q^^ {H2,iViSi + 722) is an i.i.d. AWGN 
process with zero mean and a covariance matrix propor- 
tional to the identity. Let S be the number of zeros on 
the main diagonal of Pi. Then, the matrix V2 contains 
Nr — S zero columns. Note that S = implies that no 
transmission takes place in the secondary link. In the 
sequel, we always assume that S > (which will be 
the case at low signal to noise ratio as shown in the 
simulations). 

IV. Input Covariance Matrix Optimization 

In the latter section the proposed pre-coding scheme 
does not generate any interference on the primary user 
but the transmission rate for the secondary user was not 
optimized. For this purpose, the choice of the power al- 
location of the secondary transmitter, i.e. the matrix P2, 
needs to be optimized. First, we present the most simple 
case where uniform power allocation is performed. Sec- 
ond, we introduce a power allocation which maximizes 
the individual transmission rate. In both cases we assume 
that the pre-coder has been previously adapted to satisfy 
the orthogonality conditions ([S]) or Q. 

A. Uniform Power Allocation 

For the uniform power allocation scheme the input 
covariance matrix is set to P2 = iNt and the constant 
a from Q is tuned in order to meet the condition 
Trace {V2V2) — Pmax- The rate achieved by the 
secondary user while generating zero-interference to the 
primary receiver is 



R2 = log2 



I N^ + Q-'i H 2,2V 2V^Hl2Q~ 



(12) 



B. Optimal Power Allocation 

The transmission rate for the secondary link is maxi- 
mized by adopting a power allocation matrix P2 which 
is a solution of the following optimization problem, 

^2(^2) 

(13) 



arg max 
P2 



s.t. Trace (^2^2^?) ^ 



where 



log2 



In + Q--^H2,2V2P2V^H^oQ- 



(14) 

Note that solving this optimization problem requires the 
knowledge of the covariance matrix Q, which is calcu- 
lated at the secondary receiver based on the knowledge 
of the channel H2 i- This can be done if the secondary 
receiver estimates Q and feeds it back to the secondary 
transmitter. Here, we assume a perfect knowledge of Q 
is available at the secondary transmitter. 

By definition (Eq. (fT3])). the matrix V2 is not full rank. 
Therefore, the optimization problem ([T3] ) does not have 
a simple solution. We propose a two-step optimization 
which leads to a water-filling solution. First, we define 
a new input covariance P2 such that. 



P2= {V^V2}' P2{V^V2)' . (15) 

By replacing the expression ([TSll in ([T3] ). the optimization 
problem becomes 

+ GP2G^ 



arg max 
P2 

s.t 



log2 



Trace ( P2 



(16) 



where G = Q-^H2,2V2 (^^^^2) ' G C^-^^'. The 
idea here is to solve the a priori non-trivial optimization 
problem defined by expression ( fT3l ) by introducing an 
equivalent channel matrix G to simplify the problem. 
Using G we can then apply a singular value decomposi- 
tion to the new channel such that G = EAZ^, where 
E G C^'-^^'- and Z G C^'^^' are unitary matrices, and 



the matrix A € 



contains the singular values 



r/i, . . . ,'r]^in{N^,Nt} of Under these assumptions, the 
optimal solution P2 = Z^P2Z is 



yne{l,...,Nt}, Pl{n,n) 



1 



(17) 



where p is a Lagrangian multiplier that is determined to 



N 



satisfy ^ -P2(i) J) = Pmax- Once P2 has been obtained, 

i=i _ 
then the optimal power allocation matrix lUOl is 



P2={V^V2) ZP2Z'' [V^V2) ' (18) 

The constant a in (l6]l is tuned such that the condition 
(|2]| is met for i = 2. 



V. Numerical Results 

In this section we show numerical examples to il- 
lustrate the performance of our interference alignment 
strategy. Considering the same number of antennas at 
the receiver and the transmitter, we analyze the number 
of unused singular values or free dimensions available 
for the secondary link as well as its achieved data rate. 
Recall that the primary link is interference-free, therefore 
its data rate corresponds to the single user case rate ifTOl . 

In Fig. [2] we show the number of unused singular 
values in the primary link as a function of the number 
of antennas and SNR = Note that in low SNR 

regime, the transmitter attempts to concentrate all its 
power in the best singular values leaving all the oth- 
ers unused. On the contrary, in high SNR regime the 
primary transmitter tends to spread its power among 
all its available singular values. Thus, in the first case 
the opportunistic link has plenty of free dimensions, 
while in the second one, it is effectively limited. This 
power allocation behavior has been also reported in ifTTl 
and [12|. Similarly, it is observed that increasing the 
number of antennas leads on average, to a linear scaling 
of the unused singular values. In Fig. [3l we show the 
achieved data rate of secondary link when optimal power 
allocation is implemented (i?2,optimai) as a function of 
the number of antennas and the SNR. Therein, it is 
shown that at very low and very high SNRs the data rate 
approaches zero bits/sec. In low SNR regime this effect 
is natural since detection is difficult due to the noise. 
However, at high SNRs it is due to the fact that the 
primary transmitter does not leave any unused singular 
value. Nonetheless, at intermediate SNRs, significant 
data rates are achieved by the secondary link. Note that 
increasing the number of antennas always leads to higher 
data rates for the secondary link. In Fig. IH we plot the 
data rate of the primary Ri and secondary link for the 
case of uniform i?2,uniform and optimal i?2,optimai power 
allocation. Note that in high SNR regime, the -R2,uniform 
and i?2,optimai performs similarly. It is due to the fact that 
few or even none of the singular values are left unused by 
the primary link, therefore the uniform power allocation 
does not differ from the optimal in average. In contrast, 
at intermediate SNRs the difference in performance is 
more significant for a large number of antennas. 

VI. Conclusions 

We provided a novel interference alignment scheme 
which allows an opportunistic point-to-point MIMO Unk 
to co-exist with a similar pre-existing primary link on 
the same fully utilized band without generating any ad- 
ditional interference. The proposed scheme exploits the 
fact that the licensed transmitter, while performing water- 
filling power allocation on its MIMO channel (in order to 



Primary Link (N^ = = 3) 




Figure 2. Average number of unused singular values in the primary 
link as a function of the total number of antennas and the SNR = 
^'"t'' . The SNR and the number of antennas A^r ~ Nt are assumed 
the same for the primary and secondary links 



Figure 4. Average data rate achieved by the primary Ri and 
secondary link for uniform -R2, uniform and optimal -R2, optimal power 
allocation as a function of their SNR = . The dashed lines 
correspond to Nr = Nt = 20 antennas. The solid lines correspond 
to Nr = Nt = 3 antennas. 




Figure 3. Average data rate of the secondary link when optimal 
power allocation is implemented as a function of the number of 
antennas Nr = Nt and SNR = The SNR and the number of 

antennas are assumed the same for the primary and secondary link. 



maximize its single user rate) will leave some singular 
values unused when constrained by power limitations. 
Hence, no transmission takes place along the correspond- 
ing spatial directions. We proposed a linear pre-coder 
for the opportunistic radio, which perfectly aligns the 
transmitted signal with such unused dimensions. We also 
provided a power allocation scheme which maximizes 
the data rate of the opportunistic link. Numerical results 
show that significant data rates are achieved by the 
secondary link even for a reduced number of antennas. 
Further studies will extend the novel approach to multi- 
user multi-carrier systems and the case of incomplete 
CSI. 
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